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Abstract: We describe a kinematic method which is capable of determining the over- 
all mass scale in SUSY-like events at a hadron collider with two missing (dark matter) 
particles. We focus on the kinematic topology in which a pair of identical particles is 
produced with each decaying to two leptons and an invisible particle (schematically, 
pp YY+jets followed by each Y decaying via Y iX ii'N where N is invisible). 
This topology arises in many SUSY processes such as squark and gluino production and 
decay, not to mention ti di-lepton decays. In the example where the final state leptons 
are all muons, our errors on the masses of the particles Y, X and N in the decay chain 
range from 4 GeV for 2000 events after cuts to 13 GeV for 400 events after cuts. Errors 
for mass differences are much smaller. Our ability to determine masses comes from 
considering all the kinematic information in the event, including the missing momen- 
tum, in conjunction with the quadratic constraints that arise from the Y, X and 
mass-shell conditions. Realistic missing momentum and lepton momenta uncertainties 
are included in the analysis. 
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1. Introduction 

As the Large Hadron Collider (LHC) is near completion, we will soon be able to fully 
explore TeV scale physics. Because of the naturalness problem for the Higgs boson 
in the context of the Standard Model (SM), it is strongly believed that new physics 
beyond the SM should appear at or below the TeV scale. There are many possible 
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candidates for TeV-scale physics beyond the Standard Model, giving rise to various 
experimental signatures at the LHC. If some new signal is discovered, it is vital to 
determine the masses and spins of the new particles in order to fully reconstruct the 
picture of the TeV scale. 

Some new physics will be easily identified. For example, if there is a Z' gauge 
boson accessible at the LHC, one can easily find it by looking for the resonance in the 
invariant mass distributions of its decay products, e.g., a pair of Icptons or jets. In 
general, if the decays of a new particle involve only visible particles, one can search for 
it by looking for a bump in various invariant mass combinations of the visible particles 
and the location of the bump determines the mass of the new particle. On the other 
hand, if the decays of a new particle always contain one or more invisible particles, the 
search for the new particle becomes more complicated, as there is no "bump" to look 
for. In order to detect new physics in such a case, it is necessary to understand the 
SM backgrounds very well and to then look for excesses above them. Determining the 
masses of the new particles will also be challenging since we cannot directly measure 
the energy carried away by the invisible particles. Absent good mass determinations, 
it will be difficult to reconstruct a full picture of the TeV scale even after new physics 
is discovered. 

A scenario with missing particles is highly motivated for TeV scale physics, inde- 
pendent of the hierarchy problem. If we assume that dark matter is the thermal relic 
of some weakly interacting massive particles (WIMPs) left from the Big Bang, then 
the right amount of dark matter in the universe is obtained for a WIMP mass in the 
0.1-1 TeV range under the assumption that the electrowcak sector mediates the dark 
matter — SM interaction. The dark matter particle must be electrically neutral and 
stable on cosmological time scales. If it is produced at a collider, because it is weakly 
interacting it will escape the detector without being detected, giving missing energy 
signals. In order for the dark matter particle to be stable, it is likely that there is a 
new symmetry under which the dark matter particle transforms but all SM particles 
are neutral, thereby preventing decay of the dark matter particle to SM particles. 

LEP has indirectly tested physics beyond the SM. The electroweak precision fit 
and 4-Fermi contact interaction constraints exclude new particles with masses below 
(9(TeV) if they are exchanged at tree level, unless their coupling to the SM fermions 
is suppressed. If there is a symmetry under which the new particles are odd and the 
SM particles are even, then the new particles can only contribute to the electroweak 
observables at the loop level. In this case, the bound on the mass of the new particles 
decreases by about a loop factor, m — > rajA.'K^ making the existence of new particles 
with masses of order a few hundreds of GeV compatible with the data. The message 
coming from the LEP data is that, if there is any new physics responsible for stabilizing 
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the electroweak scale, it is very likely to possess such a new symmetry. Thus, the 
cosmological evidence for dark matter together with the LEP data provide very strong 
motivation for new particles at or below the TeV scale that are pair produced rather 
than singly produced. 

Almost all the models with dark matter candidates contain additional particles 
charged under the new symmetry. At a collider, these new particles must also be pair- 
produced, and if they are heavier than the dark matter particle, they will cascade decay 
down to it. In many cases, this cascade radiates SM particles in a series of A ^ Be, 
1 — s> 2 decays, in which A and B are new physics particles while c is a SM particle. 
(In some cases, phase space restrictions force one of the new particles off-shell and 
A B*c —>■ Cdc, 1-^3 decays are relevant.) Since the final step in the chain will 
yield a dark matter particle, the typical collider signals for such a scenario will be jets 
and/or leptons plus missing energy. 

Supersymmetry (SUSY) is the most popular candidate for physics beyond the SM 
and belongs to the above category of models. In SUSY, _R-parity conservation implies 
that the Lightest Supersymmetric Particle (LSP) is stable. In most supersymmetric 
models the LSP is the lightest neutralino, which is a good dark matter candidate. 
It appears at the end of every supersymmetric particle decay chain and escapes the 
detector. All supersymmetric particles are produced in pairs, resulting in at least two 
missing particles in each event. 

Other theories of TeV-scale physics with dark matter candidates have been re- 
cently proposed. They have experimental signatures very similar to SUSY: i.e. multi- 
ple leptons and/or jets plus missing energy. For instance. Universal Extra Dimensions 
(UEDs) |T|, 1^, little Higgs theories with T-parity (LHT) [0], and warped extra dimen- 
sions with a Zs parity 0] belong to this category of models. Being able to reconstruct 
events with missing energy is thus an important first step to distinguish various sce- 
narios and establish the underlying theory. 

Of particular importance will be the determination of the absolute masses of the 
new particles, including the dark matter particle. First, these masses are needed in 
order to determine the underlying theory. For example, in the case of SUSY, accurate 
particle masses are needed to determine the SUSY model parameters, in particular 
the low-scale soft-SUSY-breaking parameters. These in turn can be evolved to the 
unification scale (under various assumptions, such as no intermediate-scale matter) to 
see if any of the proposed GUT-scale model patterns emerge. The accuracy required 
at the GUT-scale after evolution implies that low-scale masses need to be determined 
with accuracies of order a few GeV. Second, the mass of the dark matter particle, and 
the masses of any other particles with which it can coannihilate, need to be determined 
in order to be able to compute the dark matter relic density in the context of a given 
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model. Studies suggest that the required accuracy is of order a few GeV. A very 
important question is then whether or not the LHC can achieve such accuracy or will it 
be necessary to wait for threshold scan data from the ILC. One goal of this paper will 
be to find techniques for determining the dark matter particle mass at the LHC with 
an accuracy that is sufficient for a reasonably precise computation of the relic density. 

Most of the SUSY studies carried out thus far have relied on long decay chains 
of super-particles which produce many jets and/or leptons and large missing energy. 
Several kinematic variables have been proposed as estimators of the super-particle mass 
scale, such as Ht, M^ff 0, and However, these variables measure the 

mass differences between the super-particles, but not the overall mass scale. 

One possible means for determining the overall mass scale is to employ the total 
cross section. However, the total cross section is very model dependent: it depends on 
the couplings, the species of the particles being produced, e.g., fermions or bosons, as 
well as the branching fractions of the decays involved in the process. One needs to have 
already determined the spins and branching ratios for this to be reliable, a task that is 
difficult or impossible at the LHC without an ability to determine the four-momenta 
of all the particles involved in the process. To fully test a potential model, we must 
first determine the masses of the produced particles using only kinematic information. 
Once the masses are known, there are many chain decay configurations for which it 
will be possible to use these masses to determine the four-momenta of all the particles 
on an event-by-event basis. The four-momenta can then be employed in computing 
the matrix element squared for different possible spin assignments. In this way, a spin 
determination may be possible, and then the cross section information can be used to 
distinguish different models. 

In recent years there have been numerous studies of how to measure the super- 
partner masses just based on kinematic information P, ||, ^, |10], [Tl|, |T2|, |13|, 0, |16 



These studies rely on long decay chains of super-particles, usually requiring 3 or more 
visible particles in the decay chain in order to have enough invariant mass combinations 
of the visible particles. One can then examine the kinematic edges of the distributions 
of the various invariant mass combinations, and obtain the masses from the relations 
between the end points of the distributions and the masses. Many of these studies use 
the decay chain q —>■ ~^ ~^ Xi^^Q (Fig- HI) that occurs for the benchmark point 
SPSla [0, for which m^o ~ 97 GeV, 143 GeV, m^o ~ 180 GeV, ~ 570 GeV 

and rrig ~ 610 GeV, see Appendix B. The kinematic endpoints of the invariant mass 
distributions, niu, rrigu, mq^(high), and ^^^(low)^, depend on the masses of the super- 

^High and low represent the largest and the smallest values of ruqi, respectively — these masses 
are employed since it is not possible to determine the order in which the observed leptons appear in 
the chain decay. 
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particles in the decay chain through some comphcated relations 0, |T^, IT3I . If the end 
points of these distributions can be accurately determined from the experimental data, 
we can invert the relations to obtain the masses of the super-particles. 

For the decay chain of Fig. |l] and the specific model points studied, this approach 
can give a reasonable determination of the masses of the super-particles, but there is 
room for improvement. In some of the studies, it is only mass differences that are 
well determined whereas the overall mass scale is rather uncertain. For one of the 
mass points studied in |jlO|, |18| (labeled a in |jlO[), a very large number of events is 



employed and the overall mass scale uncertainty is reduced to two discrete choices, 
one corresponding to the correct solution (with rms error relative to the central value 
of order 4 GeV) and the other (somewhat less probable value) shifted by about 10 
GeV. For the mass choices labeled /3, for which the event rate is lower, there are a 
number of discrete solutions and each one has fairly large rms error for the absolute 
mass scale. However, it should be noted that in reducing the solutions to a number of 
discrete choices, not only were the locations of the kinematic edges employed, but also 
the shapes of the distributions of the mass combinations were employed. These latter 
shapes depend upon their choice of model being correct. It is possible that without 
this information there would have been a significant continuous range of overall possible 
mass scale. 

Another mass determination method is that proposed by Kawagoe, Nojiri, and 
Polesello [0. Their method relies on an even longer decay chain, ^ — > 662 — ^ X2^i^2 
£6162^2 Xi^i^2^i^2- There are five mass shell conditions and for each event there are 
the four unknowns due to the unobservable 4-momentum of the Xi- In principle, before 
introducing combinatorics and experimental resolution, one can then find a discrete 
set of solutions in the space of the 5 on-shell masses as intersections of the constraints 
coming from just five events. In practice, combinatorics and resolution complicate the 
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picture. In their actual analysis, they only fitted the gluino and the sbottom masses 
with the assumption that the masses of X2> ^^'^ Xi already known. For the 
standard SPSla point, they achieved accuracies for nig and mi of order a few GeV, 
but with central values systematically shifted (upward) by about 4 GeV. In a follow up 
study ||13[, Lester discusses a procedure for using all 5 on-shell mass constraints. For a 
relatively small number of events and without considering the combinatorics associated 
with the presence of two chains, he finds a 17% error in the determination of fn^o. 

In addition to the above studies, a series of contributions concerning mass deter- 
mination appeared in |Tl|. These latter studies focused on the SPSla point and again 
employed the kinematic edges of the various reconstructable mass distributions in the 
^ — > 662 — > X2^i^2 — > ^fci&2^2 — > Xi^i^2^i^2 decay chain to determine the underlying 
sparticle masses. Experimental resolutions for the jets, leptons and missing energy 
based on ATLAS detector simulations were employed. The resulting final errors for 
LHC/ATLAS are quoted in Table 5.1.4 of [^, assuming an integrated luminosity of 
300 fb"^ and after using both e and /i intermediate resonances (assuming they are de- 
generate in mass) . We have since verified with several ATLAS members that the quoted 
errors do indeed correspond to ilex errors ll^. The tabulated errors for /ra^^o, mj, m^o 
are all of order 5 GeV, while those for and m-g are of order 8 GeV. 

In all of the studies referenced above, the methods employed required at least three 
visible particles in the decay chain, and, in the last cases above, four visible particles 
(two 6's and two £'s). We will study the seemingly much more difficult case in which 
we make use of only the last two visible particles in each decay chain. (For example, 
the subcase of Fig. |1] in which only the X2 — > d-d-X^ portion of each decay chain is em- 
ployed.) In this case, if only the isolated chain-decays are analyzed, the one invariant 
mass combination that can be computed from the two visible 4-momenta does not con- 
tain enough information to determine the three masses {fn^o, m^and ^ri^o) involved in 
the decay chain. Thus, we pursue an alternative approach which employs both decay 
chains in the event at once. For the SPSla point, our method allows a determination 
of the masses ^^0, mj and m^o with an accuracy of ~ ±5 GeV after including both 
e and /I intermediate slepton states (again, taken to be degenerate in mass), assum- 
ing L = 300 fb~^ and adopting the ATLAS expectations for the resolutions for lepton 
momentum and missing momentum measurements. (These resolutions affect the de- 
termination of the crucial transverse momentum of the M + system. In particular, 
by looking at only the leptonic part of the decay chains we can avoid considering indi- 
vidual jet momenta, and therefore we are less sensitive to imprecise measurements for 
the individual jets.) In short, using only the leptons in the final state, we obtain an 
accuracy that is very comparable to the ~ ±5 GeV found in the LHC/ATLAS study 
referenced above for the same luminosity and very similar detector simulation. 
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In the above single-chain decay approaches, it is imphcitly assumed that appropri- 
ate cuts etc. have been employed so that both decay chains in each event involve the 
same decaying resonances, all the way back to the Tj. In our approach it is unnecessary 
to know exactly what resonances appear prior to the Xg's in the two decay chains. Thus, 
some of the X2 P^i^ events could come from direct q production and some indirectly 
from 'g production followed hj g ^ qq decay. We also do not need to tag the b quarks. 
We only need to measure to determine the transverse momentum of the X2X2 P^i^ using 
the measured lepton momenta and the measured missing momentum. Nonetheless, we 
do need to isolate a sample of events dominated by two final X2 ^ ^ ^^Xi decays. 
(Of course, it is interesting to go beyond this assumption, but we will not do so in this 
paper.) The key to mass determination using the more limited information we employ 
is to consider the whole event at once and look not for edges in masses reconstructed 
from visible momenta but for sharp transitions in the number of events consistent with 
the assumed topology after an appropriate reconstruction procedure. Further, as noted 
later, if the events we isolate do not correspond to a X2X2 P^i^ decaying in the manner 
assumed, then our procedure will yield smooth distributions in the number of recon- 
structed events, as opposed to the sharp transitions predicted if we have indeed isolated 
an enriched X2X2"P9'i^ sample with decays as presumed. 

Beginning with the general topology illustrated in Fig. ^, we employ the informa- 
tion coming from correlations between the two decay chains in the same event, and the 
missing momentum measurement. This is evident from some simple constraint count- 
ing. Each event satisfying the topology of Fig. |^ has the two invisible 4-momenta of the 
and A^'. The sum of the transverse momenta of and A^' is, however, constrained to 
equal the negative of the sum of the transverse momenta of the visible particles, leaving 
us with 6 unknowns for each event, subject to 6 on-shell mass constraints. Under the 
assumption that my = my/, mx = rux' and rri]^ = rriN', we are left with the three un- 
known masses, my, mx and rriN- Every event will be compatible with a certain region 
in the 3-dimensional {mYymxjTnx} space. Combining several events will shrink this 
region. We will show that before the inclusion of combinatorics and resolution effects 
the actual values of the masses lie at the end point of such a region. Mass determination 
after including combinatoric and resolution effects requires an examination of how the 
number of events consistent with given mass choices changes as the masses are shifted. 

In our approach, we find that it is important to not focus on the individual invariant 
mass distributions, as this would not utilize all the information contained in the data. 
Instead, we examine the events from a more global point of view and try to use all the 
kinematic information contained in the data to determine the masses of the particles 
involved.^ In the case where we have of order 2000 events available after cuts, and after 
^We note that the fact that each event defines a region in mass space was also the case for the 
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including combinatorics and resolutions for missing momentum and lepton momentum 
measurements according to ATLAS detector expectations, we achieve rms accuracies 
on my, mx and m^v of order 4 GeV, with a small systematic shift that can be easily 
corrected for. This assumes a case with significant separation between the three masses. 
This result is fairly stable when backgrounds are included so long as S/B > 2. This 
number of events and resulting error apply in particular to the SPSla point assuming 
integrated luminosity of L = 300 fb~^ and use of all i = e or Jl channels. 

The organization of the paper is as follows. In Sec. |, we give a detailed expo- 
sition regarding solving the topology of Fig. 0. In Sec. ^, we demonstrate how the 
masses of the Y, X and N particles in Fig. ^ can be very precisely determined after 
a reasonable number of events {e.g. 500) if there are no uncertainties associated with 
combinatorics or with particle and missing momentum measurement resolutions. In 
Sec. ^ we develop the very crucial strategies for dealing with the realistic situation 
where combinatorics, resolutions and backgrounds are included. Sec. ff]T| focuses on the 



effects of combinatorics and finite resolutions for the lepton and measurement missing 
momentum measurements. If of order 2000 events are present after cuts, we still find 
good accuracies for not only mass differences, but also for the absolute mass scale. 



Z ^ Y ^ X ^ N one-sided chain situation outlined earlier (except for the mass space being 4- 
dimensional) . An interesting question is whether our more general approach would determine the 
absolute mass scale in the one-sided case, as opposed to just mass differences. A detailed study is 
required. 



- 8 - 



using only the kinematic information contained in the available events. In Sec. [4.2| , we 
discuss the effects of having background events mixed with the signal events. In Sec. ||, 
we discuss two alternative scenarios: one with very different my — rrix compared to 
the first point analyzed, and one with rriAr ~ 0. In Sec. p.3| , we consider in detail the 
SPSla mSUGRA point. We summarize and present additional discussion in Sec. 0. 



2. Topology of events with missing energy 

We study the collider events with topology shown in Fig. ^ A hard hadronic collision 
produces two identical or mirrored chains. Each decay chain gives two visible particles 
and one missing particle. It will be convenient to label the 6 final outgoing particles 
from 1-6, with = 1, A^' = 2, visible particles 3 and 5 emitted from the Y chain and 
visible particles 4 and 6 emitted from the Y' chain. There are many processes which 
have this topology. For example, ti production, with t decaying to bW and W decaying 
leptonically to iu, is exactly described by this topology, so it can be studied with our 
method, except that we already know that neutrinos are (to a good approximation) 
massless. There are also many SUSY or other beyond the SM processes which can be 
described by this topology, e.g., second neutralino pair production X2X2 (through t- 
channel squark exchange) with X2 ~^ ^'^d then £ — > producing 4 visible charged 
leptons and 2 missing particles. As already noted, we require that the masses of the 
corresponding particles in the two chains be the same. They can be the same particle or 
one can be the anti-particle of the other. Or, they can even be different particles whose 
masses are believed to be approximately equal {e.g., squarks or sleptons of the first 
two generations). The visible particles do not need to be stable as long as we can see 
all their decay products and reconstruct their 4-momenta. The event can involve more 
particles (such as hadronic ISR/FSR or parent particles such as squarks and gluinos 
decaying within the gray blob on Fig. 0) as long as none of the additional particles 
lead to missing momentum. For example, the 4 leptons plus missing energy event from 
the decays of a pair of second neutralinos can be part of the longer decay chains from 
squark pair production, as occurs for the SPSla chain decay. 

It is instructive to analyze the unknowns in this topology in a more detailed manner 
than given in the introduction. In particular, we can make a distinction between 
kinematic unknowns — those in which phase space is differential — and parametric 
unknowns — Lagrangian parameters or otherwise non-kinematic unknowns on which 
the cross section has some functional dependence. For instance in the Breit-Wigner 
propagator [{(f — M^)^ + M^r^/4]~^, q is kinematic while M and F are parametric. 
Masses, including those of missing particles, are parametric unknowns (phase space 
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d^p/2E is not differential in them). Any function of an event's missing 3-momenta and 
already-known parameters is a kinematic unknown. 

Each event with the topology of Fig. ^has eight kinematic unknowns: Pn,Pn' and 
the initial state E and Pz, where we are assuming the parameters rriN and rriN' are fixed. 
Total 4-momentum conservation reduces this to four kinematic unknowns. In the nar- 
row width approximation, the mass combinations constructed from a combination of 
visible momenta and invisible momenta (which we place in the class of kinematic un- 
knowns), such as = (pi+Ps)^, are equal to the corresponding parametric unknowns, 
such as m^, and we can trade them for their corresponding parameters. Therefore, in 
the narrow width approximation, a single event is described by a volume in the six 
dimensional parameter space {my, my, mx, mx', mN,mN'}- 

If the two chains are identical or mirrors of one another and if we use the narrow 
width approximation, we can impose 3 more relations, my = my/, mx = mx' and 
mjv = mjy/, which reduces the independent unknown parameters to three. As a re- 
sult, if we know the three masses my, mx, and then (up to discrete ambiguities 
associated with multiple solutions to a quartic equation) we can solve for all the un- 
known momenta, using the measured visible momenta, and vice versa. The procedure 
is described in more detail in Appendix ^ 

If the masses are not known, we must assume values for the three masses my, mx, 
and m^v- Given a fixed M. = {my, mx, m^v} choice, for each event we obtain a quartic 
equation (for the energy of one of the invisible particles) with coefficients depending 
on the assumed masses, Ai, and visible momenta. It can have to 4 real solutions for 
the invisible energy, depending on the coefficients, and each solution fully determines 
associated 4-momenta for both invisible particles. 

Any solution with real and physically acceptable invisible 4-momenta corresponds 
to a choice for my, mx, and m^ that is consistent with that particular event. The 
points in = {my, mx, m^} parameter space that yield real solutions are not discrete; 
instead, each event defines a region in the three-dimensional mass space corresponding 
to a volume of real solutions. The region in the mass space consistent with all events, 
the 'allowed' region, will shrink as we consider more and more events. However, even 
for many events the allowed region remains three-dimensional and does not shrink to 
a point. We need to find techniques that allow us to identify the correct mass point 
given a volume in mass space consistent with a set of events. 
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3. Idealized case: Perfect resolution, no combinatorics and no 
background 



In order to understand how the mass information is contained in the kinematics, we 
start with the ideal case in which all visible momenta are assumed to be measured 
exactly and we associate each lepton with the correct chain and position in the chain 
[i.e., we neglect resolution effects and combinatorics). For illustration, we have gen- 
erated a sample of 500 events of qLQL production, with each decaying according to 
Fig. I, with Y = Y' = x^, X = X' = ilii, N = N' = Xi, and 3, 4, 5, 6 all being /I's 



of various signs. We generated our events using SHERPA pO[ versions 1.0.8 and 1.0.9 



and PYTHIA |2T|. We generated the SUSY spectrum for the mass points considered 
using SPheno 2.2.3 Details regarding the spectrum, cross sections and branching 
ratios for this point (Point I) are given in Appendix |B[ For the moment, we only need 
to note the resulting masses: 

my = 246.6 GeV , mx = 128.4 GeV , m^v = 85.3 GeV . (3.1) 

We stress however that our techniques are not specific to SUSY; we have just used 
the available tools for supersymmetric models. Thus, event rates employed are not 
necessarily those predicted in the context of some particular SUSY model. Here, we 
simply use a 500 event sample for illustration of the basic ideas. 

For simplicity, we have required all four leptons to be muons. We assume that the 
momenta of the 4 muons and the sum of the transverse momenta of the two neutralinos 
are known exactly. The only cuts we have applied on this sample are acceptance cuts 
for the muons: px > Q GeV and |?7| < 2.5. We do not consider the mass of the 
squark, therefore information from the quarks is irrelevant except that the presence 
of the quark jets typically boosts the system (Fig. ||) away from the z axis, an effect 
automatically included in our analysis. In the following, we denote a set of masses as 

= {mY,mx,rniy} and the correct set as A4a- 

Each event defines a mass region in A4 space that yields real solutions for and 
Ptv' (for which we often employ the shorter phrase 'real solutions' or simply 'solutions'). 
This region can be determined by scanning through the mass space. We then examine 
the intersection of the mass regions from multiple events. This region must contain the 
correct masses, tMa- The allowed mass region keeps shrinking when more and more 
events are included. One might hope to reach a small region near Ai^ as long as enough 
events are included. However, this is not the exemplified in Fig. ^ There, the 

3-dimensional allowed region in TVl-space is shown together with its projections on 2- 
dimensional planes. When producing Fig. ^ we discretize the mass space to 1 GeV 
grids in all three directions. As already noted, we have used the correct assignments 
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Figure 3: Mass region (in GeV) that can solve all events for the input masses 
{mY,mx,mN} = {246.6, 128.4,85.3} GeV using 500 events. 



for the locations of the muons in the decay chains. Wrong assignments will add more 
complication; this will be discussed in Sec. |[ With correct assignments, and because of 
our narrow-width and no-smearing assumptions, the correct masses A4a will result in 
at least one real p^- and pjv' solution for all events and is included in the allowed region. 
In all three 2-dimensional projections, the entire allowed region is a strip with my and 
mx close to the correct values, but left undetermined except for an upper bound. 
A lower bound is sometimes present and can be caused by the presence of events in 
which the system (Fig. 0) has a large amount of transverse momentum. The upper 
bound for m^r generally materializes using fewer events than does the lower bound. By 
examining the figures one can see that the upper bound for ttin is actually close to the 
correct uin', more generally, J^a is located near the tip of the cigar shape of acceptable 
choices in A^-space. 

An intuitive understanding of why it is that the correct mass set AiA is located at 
an end point can be garnered from Fig. ^. Any point in the mass space on the left-hand 
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Figure 4: Map between mass space and kinematic space. The nominal masses, point A, 
produces a kinematic region that coincides with the experimental region: ICa = ICexp- A 
point B inside the allowed mass region produces a larger kinematic region: /C^ D K-exp- 

side of the figure is mapped into a region of the kinematic space on the right-hand side. 
By 'kinematic space' we mean the set of observed 3-momenta of the visible particles, 
3, 4, 5, and 6. Thus, the kinematic space has much higher dimensionality than the 
mass space — the on-shell Y, X, N masses can be held fixed while changing the angles, 
magnitudes and so forth of the visible particles. Consequently, each point in mass space 
defines a volume in kinematic space. In analyzing data, the inverse mapping is to be 
envisioned. Each point in the kinematic space corresponds to a specific momentum 
configuration of the visible particles, i.e. an event. A collection of many events will 
define a region in the kinematic space. In particular, the correct set of masses, point 
A in Fig. ^, produces a kinematic region JCa that coincides with the experimental one, 
K-A = JCexp, as long as the number of experimental events is large enough so that all 
the allowed region is populated. Any shift away from A will generally not allow one or 
more kinematical observables associated with the visible particles to occupy a region 
close to the boundary of JCexp] i-e. such a shift will generally exclude a region of the 
actually observed kinematical space. 

A mass point other than J^a produces a region different from JCexp- If it does not 
cover the entire JCexp, this means that some events will not have yielded real pn and 
Pn' solutions and, therefore, the mass point does not appear in the final allowed mass 
region. On the other hand, there can be mass points which produce larger kinematic 
regions encompassing the entire JCexp region. These mass points yield real solutions for 
all events and hence belong to the final allowed region. This kind of point is exemplified 
by point B in Fig. ^. If we shift such a point in the mass space by a small amount, 
M.B ^ M.' = J^B + <5A^, the resulting kinematic region still covers JCexp- In this case. 
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Ai' still yields real solutions for all events. Thus, point B, which produces a region 
larger than JCexp, has the freedom to move in many directions because it lives inside the 
allowed region rather than on its boundary. On the other hand, the correct mass point 
A, which produces exactly ICexp, has the least freedom to move. In short, locating the 
correct mass point Ma can be viewed as a kind of generalization of the 'edge' method 
which employs sharp edges in certain invariant mass combinations constructed from 
the visible momenta. Our method matches the whole boundary of the allowed region 
in the high-dimensional kinematic space of the visible momenta. 

Of course, using the "tip" of the allowed mass region is not applicable in the 
realistic case where experimental resolutions and combinatorics are included, not to 
mention the possible presence of background events. In particular, some of the events 
generated after including these effects will be inconsistent {i.e. not yield real solutions 
for and pn>) with the correct mass set Ma aiid so this point will not be contained in 
the Ad volume obtained if all events are considered. We must find more sophisticated 
methods to identify the correct mass point. Nevertheless, understanding the idealized 
case provides useful guidance for understanding how to deal with the more complicated 
realistic situations. 



4. Inclusion of combinatorics, finite resolutions and backgrounds 

In this section we discuss the more realistic case with finite resolutions, combinatorics 
and backgrounds. We first discuss the effects from finite resolutions and combinatorics 
and later we will include the backgrounds. For the moment, we continue to employ 
the spectrum associated with the SUSY Point I, as specified in Appendix 0, with 
{my,mx,mjv} = {246.6,128.4,85.3} GeV. 

4.1 Finite resolution effects and combinatorics 

Experimental effects related to smearing and combinatorics will deform or even kill 
the allowed mass region. In particular, since the correct mass point is located at the 
endpoint, it is most vulnerable to any mismeasurement. This can be seen in Fig. |^, 
which corresponds to 500 events for the same mass point as Fig. |^. The difference is 
that we have: i) added smearing; ii) considered all possible combinatoric assignments 
for the location of the muons in the two decay chains; and iii) included the finite widths 
of the Y and X intermediate resonances. We smear muon momenta and missing pt 
using the low-luminosity options of the ATLAS fast simulation package ATLFAST as 
described in Sees. 2.4 and 2.7 of pSf. Very roughly, this corresponds to approximately 



Gaussian smearing of muon momentum with width ~ 3%/pt and of each component 
of missing momentum p™^^ with width ~ 5.8 GeV. We note that we are not directly 
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sensitive to the resolution associated with individual jet momentum measurements; 
uncertainties in the determination of individual jet momenta are, of course, reflected in 
the uncertainty of the determination of p^^^^ as we shall shortly review. Our approach 
is only sensitive to p™'^** uncertainties because we do not look at the jets associated with 
the chain decays prior to arriving at the X2X2 pair. We only need the net transverse 
momentum of the X2X2 P^i^ a whole, and we determine this in our procedure as 
YlePT +Pt^'^'^- Thus, in our analysis the errors from smearing derive entirely from the 
uncertainties in the lepton and missing momentum measurements. The fact that we 
don't need to look at individual jets is, we believe, and important advantage of our 
approach to determining the ^ ^ind Xi masses. Of course, once these masses have 
been determined, the edge techniques, which fix mass differences very accurately, can 
be used to extract the g and q masses. 

We summarize the missing energy procedure as described in Sec. 2.7 of [23] in a 
bit more detail. The missing transverse energy EJ^^^^ is calculated by summing the 
transverse momenta of identified isolated photons, electrons and muons, of jets, 6-jets 
and c-jets, of clusters not accepted as jets and of non-isolated muons not added to 
any jet cluster. Finally, the transverse energies deposited in cells not used for cluster 
reconstruction are also included in the total sum. Transverse energies deposited in 
unused cells are smeared with the same energy resolution function as for jets. From the 
calculation of the total sum E^'^ the missing transverse energy is obtained, i?™'^^ = E^^ 
as well as the missing transverse momentum components, p™'*^*^ = —p'^'^ and p'^^'^'^ = 

_ obs 

■ 

For combinatorics, we assume no charge misidentification. Then, there are 8 in- 
dependent possible combinatoric locations for one event, which can be reduced if one 
muon pair is replaced by an electron pair. If any one of these 8 possibilities yields a 
real solution (after including smearing/resolution as described above), we include the 
A4 point in our accepted mass region. 

As regards the resonance widths, these have been computed within the context of 
the models we have considered, as detailed in Appendix p. In our Monte Carlo, the 
mass of a given X2 ^ resonance is generated according to a Breit Wigner form using 
the computed width. Although there will be some model dependence of the widths 
in that they might differ between the SUSY models employed as compared to a little- 
Higgs model, the widths for these weakly interacting particles are all much smaller than 
detector resolutions in both models {e.g. of order a few hundred MeV in the SUSY 
models). This is again an advantage of our approach since we never need to know where 
on the Breit- Wigner mass distribution of the ^ and q resonances a given event occurs. 
We only need the net transverse momentum of the X2X2 system as determined from 
YliPr + Pt'^*^- Also, for the moment we will focus on events with four /x's in the final 
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state and so both sleptons in the two decay chains will be /x's. When we come to the 
SPSla mSUGRA point, we will discuss combining results for the Jljl, ee and Jle decay 
channels. Even in this case, we analyze final states with definite lepton composition 
(4/i, 4e or 2e2/i) separately and do not need to worry about whether the Jl is closely 
degenerate with the e (although it in fact is). If there is significant non-degeneracy, 
that would emerge from our analysis. However, to get final errors on the i mass as low 
as ~ 5 GeV, degeneracy of the Jl and e must be assumed (and of course is predicted 
in the model). If in some other model, the Jl and e are not degenerate, then errors on 
these individual masses will be of order ~ 10 — 12 GeV, but errors on m^o and m^o 
will only be slightly larger than the ~ 5 GeV level since the different channel analyzes 
can be combined for their determination. 

The effects of both wrong combinations and smearing are manifest in Fig. |^: wrong 
combinations increase the chance that a given event can be 'solved' ^ and therefore 
broaden the allowed region for low itln. On the other hand, the allowed region has 
shrunk in the rriN direction with the new upper bound corresponding to a much smaller 
value. This can also be understood by using Fig. ^: some events near the boundary of 
)Ca can be resolution-smeared to a location outside of )Ca, which renders JCexp larger 
than ]Ca- Thus the correct mass point A is removed from the allowed mass region. 
For point B which corresponds to a larger kinematic region, if the fiuctuation is small 
enough, JC^xp is still within Kb and therefore does not disappear. Of course, if the 
smearing is large, the entire allowed region can be eliminated. The effect from back- 
ground events, as considered in the next subsection, will be similar. Since background 
events are produced by some completely different processes there is no reason to expect 
that multiple background events can be solved by the assumed topology with a given 
choice of M.. Thus, background events tend to reduce the allowed region. 

From the above observation, one concludes that allowed mass region in general 
does not exist and, even if it exists, we can not read directly from it the correct masses. 
Some other strategy must be employed. An obvious choice is to examine the number of 
solvable events for various given masses. We can not simply maximize the number of 
solvable events and take the corresponding masses as our estimate — such a procedure 
would still favor low m^v values. Instead, we choose to look for the mass location where 
the number of solvable events changes drastically. This kind of location is most easily 
illustrated in one dimension. For example, in Fig. ^ a, we fix my and mx to the correct 
(input) values, and count the number of solvable events as a function of m^. (In this 
figure and the following discussion, we use bin size of 0.1 GeV). A sudden drop around 
the correct ttin is obvious. Similarly, in Figs. ^ b and ^ c we have fixed my and rriN 

•^We define a 'solved' event to be an event such that the given {my, uixjItin} choices yield at least 
one solution to the final quartic equation that leads to physically allowed values for pN and pn'- 
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Figure 5: The allowed mass region (in GeV) with smearing and wrong combinatorics. 



(a) fit (b) fit (c) fit 
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Figure 6: One-dimensional fits by fixing the other two masses at the correct values. 



{mx and rriN) and also see clear "turning points" near the correct mx {fny) mass. To 
pin down where the turning points are located, we fit Figs. ^ a and || c to two straight 
line segments and take the intersection point as the turning point. 
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We can not fix a priori two of tfie masses to tfie correct values since tliey are 
unknown. On tlie otlier liand, to searcli for tfie sharpest turning point directiy in the 
3-dimensionai space is numericaffy non-trivial. This observation motivates us to obtain 
the masses from a series of one-dimensional fits. We start from some random set of 
masses and carry out a recursive series of one- dimensional fits to the number of solved 
events as a function of m^r, mx or my holding {my,mx}, {my,mAr}, or {mx,mN} 
fixed, respectively. Each such one dimensional fit gives us a sharp turning point that 
is used to set an updated value for m^, mx or my, respectively. We use this new 
value in performing a fit for the next mass in the sequence in the next step. One might 
hope that this procedure will converge to the correct mass values, but in practice, even 
though the procedure passes through the correct mass point, the fitted masses keep 
increasing and the recursion does not stabilize at the correct mass point. However, as 
we will see, there is a simple way to get the masses out of the fits. 

Having discussed the main ingredients of the method, we present a specific proce- 
dure for obtaining the masses. The procedure is applied to a data sample corresponding 
to 90 fb"^ at the LHC, using the event rates and branching ratios obtained for the SUSY 
Point I as detailed in Appendix which, in particular, gives the same masses as those 
employed in Sec. ^. {my, mx, w^at} = {246.6,128.4,85.3} GeV. Taking into consider- 
ation the decay branching ratios, the number of events is roughly 2900. In order to 
mimic reality as much as possible, experimental resolutions and wrong combinatorics 
are included. To reduce the SM background, we require that all muons are isolated 
and pass the kinematic cuts: 



With these cuts, the four-muon SM background is negligible. The number of signal 
events is reduced from 2900 to about 1900. 

The procedure comprises the following steps: 

1. Randomly select masses my > mx > that are below the correct masses (for 
example, the current experimental limits). 

2. Plot the number of solved events, Nevt-, as a function of one of the 3 masses in the 
recursive order mAr, mx-, rny with the other two masses fixed. In the case of my 
and m AT, we fit Ng^t for the plot with two straight lines and adopt the mass value 
at the intersection point as the updated mass. In the case of mx, the updated 
mass is taken to be the mass at the peak of the N^yt plot. 

A few intermediate one-dimensional fits are shown in Fig. ^ 
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3. Each time after a fit to rriN, record the number of events at the intersection 
(sometimes called the turning point) of the two straight lines, as exemplified in 
Fig. ^ a. This event number at the turning point will in general be non-integer. 

4. Repeat steps ^ and |[ The number of events recorded in step |^ will in general 
increase at the beginning and then decrease after some steps, as seen in Fig. p. 
Halt the recursive procedure when the number of (fitted) events has sufficiently 
passed the maximum position. 

5. Fit Fig. H to a (quartic) polynomial and take the position where the polynomial 
is maximum as the estimated rriN. 

6. Keep m^v fixed at the value in step |^ and do a few one-dimensional fits for my 
and mx until they are stabilized. Take the final values as the estimates for my 
and mx- 



(a) fit (b) fit (c) fit 




m„ (GeV) mx (GeV) niy (GeV) 



Figure 7: A few steps showing the migration of the one dimensional fits. The middle curve 
in each plot corresponds to masses close to the correct values. 

A deeper understanding of our procedure can be gained by examining the graphical 
representation of the steps taken in the (my,mAr) plane shown in Fig. |^. There, we 
display contours of the number of (fitted) events after maximizing over possible mx 
choices. The contours are plotted at intervals of 75 events, beginning with a maximum 
value of 1975 events. As we go from 1975 to 1900 and then to 1825 events, we see 
that the separation between the contours decreases sharply and that there is a 'cliff' 
of falloff in the number of solved events beyond about 1825 events. It is the location 
where this cliff is steepest that is close to the input masses, which are indicated by 
the (red) star. The mass obtained by our recursive fitting procedure is indicated by 
the (blue) cross. It is quite close to the actual steepest descent location. It is possible 
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that use of the contour plot by visually picking the point of steepest descent might 
also yield an accurate mass determination comparable to or possibly even superior 
to that obtained (and specified in detail below) using the recursive fitting procedure. 
Roughly, the steepest descent point corresponds to the point where the magnitude 
of in mass space is maximized. Unfortunately, even after some smoothing, the 
second derivative is quite 'noisy' and therefore not particularly useful in a local sense. 
The one-dimensional fits give us a quick and intuitive way to find this maximum, and 
the associated recursive procedure has the advantage of being insensitive to statistical 
fluctuations in the number of events at a single point. Of course, if one has the computer 
power, probably the most accurate procedure would be to directly fit the 3-d Nf^^t vs. 
{my, mx, rnisi} histogram. Fig. || is constructed from a 1-d projection of the 3-d space, 
and has therefore lost some information. 

Following the recursive fitting procedure, the final values for the masses are deter- 
mined to be {252.2, 130.4, 85.0} GeV, which are all within a few GeV of the actual 
input values, {246.6, 128.4, 85.3} GeV. The procedure is empirical in the sense that 
many of the steps could be modified and improved. In particular, above we adopted 
the criterion that the correct masses maximize the number of events at the turning 
points in the rriN fits, which is justified by Fig. |^ a. Instead, we might opt to maximize 
the number of events in the mx fits shown in Fig. |^ b. One could also change the order 
of fits in step ^ and change the fit function from straight lines to more complicated 
functions, etc. We have tried several different strategies and they yield similar results. 
Finally, one could simulate the signal for a mass point and directly generate Fig. |^, 
changing the masses until we get the best possible fit to the data; but, this is very 
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Figure 9: Contours for the number of solved events in the rriAr ~ my plane with 2000 
events. The number of events is the maximum value obtained after varying mx- Contours 
are plotted at intervals of 75 events, beginning with a maximum value of 1975. The red star 
is the position for the correct masses and the blue cross is the position of the fitted masses. 
The green dots correspond to a set of one-dimensional fits. 

computationally intensive. 

The recursive procedure does not provide an easy way to evaluate the errors in the 
mass determination. For this purpose, we generate 10 different data samples and apply 
the procedure to each sample. As above, each sample corresponds to 1900 experimental 
data points after cuts. Then, we estimate the errors of our method by examining the 
statistical variations of the 10 samples. This yields central masses and rms errors of 

my = 252.2 ± 4.3 GeV, = 130.4 ± 4.3 GeV, = 86.2 ± 4.3 GeV. (4.2) 
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The statistical variations for the mass differences are much smaller: 

my -mx = 119.8 ± 1.0 GeV, - m^v = 46.4 ± 0.7 GeV. (4.3) 

Compared with the correct values, A4a = {246.6, 128.4, 85.3} GeV, we observe small 
biases in the mass determination, especially for the mass differences, which means that 
our method has some "systematic errors". (The biases will, of course, depend upon the 
particular functions employed for the one dimensional fits — our choice of using straight 
lines is just the simplest.) One technique for determining the biases is to perform our 
analysis using Monte Carlo data. In particular, one could examine the plots of number 
of 'solved' events vs. test mass as obtained from the data vs. those obtained from a 
Monte Carlo in which definite input masses (which are distinct from the test masses 
employed during our recursive procedure) are kept fixed. One would then search for 
those input masses for which the distributions of the solved event numbers from the 
Monte Carlo match those from the data. Knowing the underlying Monte Carlo masses 
as compared to the extracted masses would allow us to subtract the differences, thereby 
removing the biases. This procedure would not appreciably change the errors quoted 
above. We beheve that the biases are mainly a function of the underlying masses and 
broad kinematic event features. However, there may be some weak dependence of the 
biases on the actual model being employed. Within the context of a given, e.g. SUSY, 
model, the bias can be quite accurately determined. 

In the above error estimation, we have neglected the uncertainties coming from 
varying the choice of the starting point in mass space used to initiate the recursive 
sequence of fits. This may introduce an error for the absolute mass scale of order the 
step size around the correct masses. For the masses chosen, it is about 1 GeV and 
much smaller than the uncertainties from varying data samples. 

The reader may be surprised at the small size of the errors quoted above given that 
the error in the measurement of the missing momentum of any one event is typically of 
order 5 GeV or larger. The explanation is similar to that associated with understanding 
the small errors for the edge locations in the edge approach. In the edge approach, the 
location of the edge for some mass variable myis is obtained by fitting data obtained at 
several mvis values. Each such data point has many contributing events and the average 
value will obviously have much smaller error than the value for any one contributing 
event. The fit to the edge will further reduce sensitivity to individual events. In our ap- 
proach, the edge in the distribution of Ng^t as a function of one of the trial masses (mjv, 
mx or my) will similarly be an average over many events and the uncertainty of the 
location of this edge will be much smaller than the uncertainties in the measurements 
of the lepton momenta and missing momentum of any one event. 
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Figure 10: Fits with 1900 signal events (after cuts) and an equal number of background 
events. Separate numbers of signal (blue) and background (red) events are also shown. 



4.2 Backgrounds 

For the point we have chosen with a 4 muon + missing energy final state, the background 
is negligible. We examined backgrounds arising from ZZZ, ZWW, tt, ttZ, ttbb, and 
bbbb. Muons from bottom and charm decays are never very hard nor isolated, and can 
be easily separated from the signal with basic isolation criteria. Tri-boson production 
simply has tiny cross sections, especially after requiring all-leptonic decays. 

Thus, we must 'artificially' introduce background in order to see what its effect 
might be on our procedures. For this purpose, we generate ti events, where the VF's 
decay to muons. We require that the b quarks decay to muons, but do not require them 
to be isolated. In many ways, this is a near-worst case background since it has a similar 
topology aside from the final 6 — > /i + . . . decays. However, the missing neutrinos imply 
that the missing momentum may be significantly different. As noted, this is not a 
realistic background as it could be removed by simple isolation cuts on the muons. 

Adding a number of background events equal to the number of signal events, i. e. 
1900 events after cuts, we repeat the one- dimensional fits. A typical cycle around the 



correct masses is shown in Fig. For comparison the numbers of solvable signal events 
and background events are also shown separately. The effect of background events is 
clear: the curve for solvable background events is much smoother around the turning 
point, and therefore smears but does not destroy the turning point. Although we are 
considering one specific background process, this effect should be generic, unless the 
backgrounds happen to have non-trivial features around the turning points. Neverthe- 
less, due to the fact that there are 8 possible combinatoric muon locations, the chance 
that a background event gets solutions is quite large and they do affect the errors and 
biases of the mass determination. This can be seen in Fig. O, in which we have used 
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Figure 11: determination with different background-signal ratio. The dashed horizontal 
line corresponds to the correct niN- 



the same 10 sets of signal events as in the previous subsection, but varied the number of 
background events according to the ratio B background) / S{ignal) = 0,0.1,0.2,0.5,1. 
We observe increases in both the biases and variations about the central values. For 
B/S > 1, the maximum in the rriN determination is obscured or even lost and we start 
to get random results. For B/S < 0.2, we are close to the B = results. 

It is important to emphasize that the above analysis is pessimistic in that it assumes 
that we do not understand the nature/source of the background events. One procedure 
that could significantly improve on uncertainties associated with the background would 
be to Monte Carlo the background or use extrapolations of measured backgrounds [e.g. 
those as measured before the cuts that make the signal a dominant as compared to a 
small component of the observed events) and then apply our recursive procedure to 
the known background and at each stage subtract off the background contribution to 
a given plot of the number of events vs. m^r, my or rrix- After such subtraction, 
the recursive procedure will yield essentially identical results to that obtained in the 
absence of background unless the background itself is not smooth in the vicinity of the 
'turning' points. 

The importance of finding cuts that both select a given topology and minimize 
background is clear. If it should happen that the we assume the wrong topology for 
the events retained, then our analysis itself is likely to make this clear. Indeed, events 
with the "wrong" topology would almost certainly yield a smooth distribution in plots 
of retained event number vs. any one of the masses of the resonances envisioned as 
part of the wrong topology. It is only when the correct topology is employed that sharp 
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steps will be apparent in all the event number vs. resonance mass plots. 

Another important situation to consider is that in which it is impossible to find 
a set of cuts that isolates just one type of decay topology, so that there are several 
signal processes contributing after a given set of cuts. However, it is quite easy to find 
situations where there are different signal processes yielding very similar final decay 
topologies, all of which would be passing through our analysis. One must then look for 
additional tricks in order to isolate the events of interest. In some cases, this is possible 
on a statistical, but not event-by-event basis. The SPSla SUSY point provides an 
interesting example that we will consider shortly. 



5. Other processes and mass points 



Our method is generic for the topology in Fig. 0, and in particular is not restricted 
to the SUSY process we have considered so far. The statistical variations and biases 
probably do depend to some extent on the process. For example, if the visible particles 
5 (6) and 3 (4) are of different species, the number of wrong combinatorics will be 
reduced and we would expect a better determination of the masses. On the other 
hand, if one or more of the visible particles are jets, the experimental resolution and 
therefore the statistical error will be worse than in the 4-lepton case. 



5.1 Changing relative mass differences in the chain 

The errors in the mass determination also depend on the mass point, especially the 
two mass differences, Amyx = my — rrix and AmxN = i^x — ttin. In Fig. [1^, a set 
of one-dimensional fits are shown for mass point Ai = {180.8, 147.1, 85.2} GeV (which 
we label as Point II). We will assume 2000 events after cuts, very similar to the 1900 
events remaining after cuts for Point I. Point II differs from Point I in that for Point II 
AniYx < AnixN, while for Point I Amyx > Auixn- The double peak structure in the 



Point II nix fit (Fig. 12 b) is evident. The curve to the right of the turning point in 
Fig. |T2| c is also "bumpy" compared with Fig. |^ c. These features are induced by wrong 
combinatorics. In the process we consider, all visible particles are muons, so they could 
be misidentified as one another and still yield solutions. Roughly speaking, Amyx and 
AmxN determine the momentum of the particles 5 (6) and 3 (4) in Fig. 0, respectively. 
Therefore, the chance that a wrong combinatoric yields solutions is enhanced when, for 
example, Amyx is close to the correct value of Am^Ar. When the two mass differences 
are close to each other, the turning point is smeared. Nonetheless, with 2000 events 
after cuts, the errors obtained for the masses are similar to those obtained for Point I. 
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Figure 12: One-dimensional fits for mass point {180.8, 147.1, 85.2} GeV. 



5.2 Small LSP mass 

Another interesting case is that of m^v being zero or very small. As for the previous 
case, we have arbitrarily used a sample of 2000 events after cuts. Because the one- 
dimensional fits proceed in the direction of increasing masses, we will miss the correct 
masses even when we start from m^v = 0. Since we always fit the ttin plot to two 
line segments, it will never yield rriN = 0. However, we can distinguish this case by 
looking at the peak number of events in the mx fits. For example, considering mass 
point {199.4, 100.0, 0.1} GeV (which we call Point 111), we start from mx = 80.0 and 
mjv = 0.0 and fit the masses in the order: my — > mx — > rrix- The first few fits yield 

{205.0, 80.0, 0} ^ {205.0, 101.5, 0} ^ {205.0, 101.5, 24.6} ^ ■ • • 

After only two steps, the Y and X masses are adjusted close to the correct values. 
Examining the peak number of events in the mx fits (Fig. |13D, we find that the number 
is maximized in the first mx fit. This is clearly different from previous cases where the 
number of events always increases for the first few mx fits (see Fig. ^b), and indicates 
that rTT-AT is near zero. 

5.3 The SPSla Point 

It is desirable to compare directly to the results obtained by others for the SPSla SUSY 
parameter point. We perform the analysis using the same AfiXiXi final state that we 
have been considering. For the usual SPSla mSUGRA inputs (see Appendix Q the 
masses for Y = xt, X = JIr and N = Xi (from ISAJET 7.75) are 180.3 GeV, 142.5 GeV 
and 97.4 GeV, respectively. This is a more difficult case than Point 1 considered earlier 
due to the fact that the dominant decay of the X2 is X2 ~^ '^^i- The branching ratio 
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Figure 13: Peak number of events in mx fits for mass point {199.4, 100.0, 0.1} GeV. 



for X2 ~^ IW'R is such as to leave only about 1200 events in the 4/iXiXi ^"^sl state after 
L = 300 fb~^ of accumulated luminosity. Cuts reduce the number of events further 
to about 425. This is too few for our technique to be as successful as for the earlier 
considered cases. After including combinatorics and resolution we obtain: 

my = 188 ± 12 GeV, = 151 ± 14 GeV , m^v = 100 ± 13 GeV . (5.1) 

In Fig. n, we give an SPSla plot analogous to Fig. ^ Errors are determined by 
generating many such plots for different samples of 425 events. Note the vertical scale. 
The change in the number of events as one varies m^v is quite small for small event 
samples and this is what leads to the larger errors in this case. 

In principle, we must also take into account the fact that the xH ~^ 'T'^i decays 
provide a background to the purely muonic final state. The dominant decay X2 
has a branching ratio that is a factor of ~ 14 times larger than that for xH ^ ^^R 
The Ti will then decay to tx\- If both r's then decay to /iz/z7, then xH ~^ ''"^1 events 
will be likely to contaminate the xH ~^ f^J^R sample. Fortunately, this contamination 
is not huge. The relevant effective branching ratios for X2X2 ~^ ttitti '^^^x\x\ ^ind 
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0.85, 



(5.3) 



"'This is, of course, due to the fact that X2 prefers to couple to left-handed slepton components, 
which are significant for the ri . 
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Figure 14: Fitted number of events at the turning point as a function of ttt-at for the fits 
for the SPSla case. 



respectively. The contamination levels from these backgrounds are further reduced by 
factors of ~ 5 for the X2X2 ~^ ttitti final state and by ~ 2 for the X2X2 tti^Jxr final 
state after imposing the simple cuts of Eq. (|4.1|) (due to the softer nature of the /x's 
coming from the r decays), implying contamination at about the 3.6% and 40% levels, 
respectively. Clearly, it is important to reduce this level of contamination given that 
m^r^ is smaller than fnj^ by about 15 GeV and so, to the extent that events containing 
X2 TTi decays remain in our sample, they might contribute additional structures 
to our plots of number of events vs. mass. This reduction can be accomplished on a 
statistical basis using a further trick analogous to that discussed (but not, we believe. 



actually employed) in Ref. |TD]. They note that the decay sequences X2 ^ ^^Xi ^"^^ 



fi e~Xi are unique to X2 ^ Thus, when considering just the one-sided decay 
chain situation one can subtract off (on a statistical basis, i.e. after many events) the 
X2 background by 



ml 



Nixl 



/i/iX?) - Nixl 



/^ex?) 



(5.4) 



where N is the number of 'solved' events as a function of one of the unknown on-shell 
masses. In our case, where both chain decays are considered simultaneously, we have 
4/iXiXi states arising from X2X2 — > t^TiT^Ti" decays and X2X2 ~^ ''"^^i^/^^/^ij decays 
in addition those from our X2X2 ~^ I^^V^'rI^'^V^ signal. To subtract off the background 
SUSY events from the former two decay chains, we can employ the following subtraction 
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(where the initial X2X2 ^^'^ final XiXi implicit) 

Niii^Jilli^Jil ^ 11+11-11+11-) = N{ii+ii-ii+ii-) - N{e+ii-ii+ii-) + N{e+e+ii-ii-) 

— N{fj,+ fj,^ fj,+ fj,^) — - N{e+ fi^ fi+fi^) + N{e- fi+ fi^ fi+) 

+A^(//+e"e+e") + A^(/x"e+e"e+) 

1 
'2 



+^ [iV(e+e+/i"/i") + N{e-e-fi+fi+ 



(5.5) 

where the latter form is likely to have smaller statistical error. An experimental indi- 
cator of the sensitivity to statistics could be gained by examining the different possible 
equivalent subtractions, of which only two are indicated above. If one were happy to 
ignore the 3.6% contamination from t+t^t+t+ decays one could then use 

a simpler form to subtract off the dominant contamination from X2X2 ~^ '^^^ t^^^V^ 
decays, namely 

N{ii+il+ii+'jl+ 11+11- 11+11-) ~ Ndi+ii'ii+ii-) - N{e+ii-ii+ii-) 

~ N{^+li-^+^-) - i [N{e+^-^+^-) + N{e-^+^-^i+)] . 

(5.6) 

We have not actually performed this kind of analysis using any of the possible subtrac- 
tions to see how well we do, but we expect that the net background contamination will 
be equivalent to B/S < 0.1, a level for which our techniques work very well and the 
errors quoted earlier for the SPSla point using the A/i final state will not be increased 
by very much. 

Of course, the same analysis as performed for the 4/i final state can also be used for 
the 2/i2e and 4e final states. Combinatorics are less of an issue for the 2//2e final state 
than for the A/i and 4e final states. The 4e-channel event number is essentially the same 
as the 4//-channel event number and the 2//2e-channel event number is roughly twice 
as large. Combining all channels (as appropriate if the e has mass very close to the 
/i, as predicted by the model), one obtains a total of about 1700 events and ~ 5 GeV 
errors for our m^o, mj and m^o determinations. 

Of course, as the observant reader may have noticed, to get 1700 events requires 
running at high luminosity, whereas the simulations referenced so far have employed the 
p™'^® resolution expected at low- luminosity running. The p™'^^ low-luminosity resolution 
is about 5.8 GeV and that at high luminosity is about 11.6 GeV. However, we have 
argued that perhaps one is not all that sensitive to this resolution when considering 
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a large collection of events and looking for the location of an edge in the number of 
reconstructed events. We have used the SPSla point to test this hypothesis by repeating 
our analysis using high-luminosity running resolutions. 

The results confirm our hypothesis. First, the worse resolution results in our ac- 
cepting somewhat more events than previously, roughly 480 events (this is the average 
for the 10 Monte Carlo "experiments" employed) in the 4/i channel for the same cuts 
(and L = 300 fb~^). The resulting mass determinations obtained using the 10 inde- 
pendent Monte Carlo experiments are 

my = 187 ± 10 GeV , = 151 ± 10 GeV , mAr = 98 ± 9 GeV, (5.7) 

where the errors are, as always, rms errors. In short, we get even smaller errors than 
for low-luminosity running. After combining the 4yU, 4e and 2/x2e channels assuming 
e~]l degeneracy our mass determination errors are slightly above 4 GeV. 

6. Summary and Discussion 

For any theory that simultaneously provides a solution of the hierarchy problem and a 
dark matter particle as a result of a symmetry guaranteeing its stability, implying pair 
production of its heavier partners, the relevant LHC events will be ones in which the 
heavier partners are pair produced, with each chain decaying down to largely visible 
SM particles and the dark matter particle, which we denote by N. In many interesting 
cases, towards the end of each such chain 2 visible SM particles emerge along with 
the invisible dark matter particle, e.g. Y fiX fifiN, with the preceding parts 
of the decay chains giving rise to jets. In other cases, two Y particles are directly 
produced and initiate 2 such chain decays. In this paper, we have developed a highly 
effective technique for using the kinematic information in a typical event containing two 
Y — > /iX — ^ /i/iiV decay chains to determine not just the mass differences in the chain 
decay, but also the absolute mass scale, using only the measured /i momenta and overall 
visible and missing transverse momenta. Since we use purely kinematic information, our 
mass determination does not require any assumptions regarding particle spins, shapes 
of distributions, cross section and so forth. Further, our procedure works whether or 
not we know the topology of each of the chains that precedes the Y ^X fifiN 
stage. This can be a big advantage. For example, in the supersymmetry context this 
allows us to combine g and q initiated chains. 

In our study, we have included resolution smearing for muon momenta and miss- 
ing momentum as incorporated in the ATLFAST simulation program. We have also 
included full combinatorics appropriate to the jets + A/iNN final state. Assuming of 
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order 2000 events after cuts and ATLFAST resolutions appropriate to low-luminosity 
running, we have found statistical errors of order 4 GeV for the individual Y, X and 
N masses, assuming a reasonable background to signal ratio, B/S < 0.5. There is also 
a small systematic bias in the masses extracted. However, this bias can be removed 
using Monte Carlo simulations once the masses are fairly well known. The appropriate 
procedure is described in Sec. [4.1| . We have not yet performed the associated highly 
computer intensive procedure, but believe that the systematic biases can be reduced 
below 1 GeV (a residual that we think might arise from possible model dependence of 
the kinematic distributions). 

As a particular point of comparison with the many earlier studies that use the mass- 
edge technique, we have examined the standard SPSla point. Following our procedure 
we are left with about 1920 events (averaging over 10 Monte Carlo "experiments") in 
the jets + 4fi, jets + 2e + 2fi and jets + 4e channels after cuts assuming an integrated 
luminosity of 300 fb^^ and employing resolutions appropriate to high-luminosity run- 
ning. The errors on m^, m^and m^o are all between 4 GeV and 5 GeV if Jl and e mass 
degeneracy is assumed. The previous mass-edge studies make this same assumption 
and employ all the final SM particles of the full g ^ bb ^ bbx% bbii bbiixi decay 
chain but examine only one chain at a time. Only one of these mass-edge studies claims 
an accuracy (~ ±5 GeV for m^o, and ttl^o) for the same channels and integrated 
luminosity that is competitive with the small error we obtain. 

By comparing the SPSla results obtained for high-luminosity resolutions to those 
for this same point using low- luminosity resolutions (as summarized in the previous 
section) we found the important result that the accuracy of our mass determinations was 
very little influenced by whether or not we employed low- or high-luminosity resolution 
for p^^^^, the latter being essentially twice the former. That our ability to locate the 
"edge" in a plot of the number of reconstructed events, N^vt, as a function of the 
test value of, say, m^o, is not noticeably affected by a factor of two deterioration in 
resolution for p^^^^ is a sign of the robustness of our approach. 

Accuracies of order 4 — 5 GeV for the masses of new-physics particles will yield 
an accurate determination of the TeV-scale underlying parameters of the associated 
new physics model. The latter accuracy will, in turn, typically yield reasonably precise 
evolved values for the model parameters at any higher scale {e.g. the coupling constant 
unification scale in SUSY) where they might follow a meaningful pattern that would 
determine the more fundamental structure of the new physics theory. Further, an 
accuracy of order 4 — 5 GeV for the dark matter particle mass will in many cases allow 
a sufficiently accurate calculation for the dark matter density from annihilations in the 
early universe as to allow a meaningful comparison with the very accurate observed 
value. In some cases, the dark matter particle coannihilates with another particle of 
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only slightly larger mass. We will be exploring the extent to which the mass of the 
coannihilation partner could be determined in such a situation. For the moment, we 
can only claim the 4 — 5 GeV kind of error on individual masses when mass differences 
are reasonably substantial (and the number of events after cuts is of order 1700 to 
2000). 

A 'fun' example that we hope our experimental colleagues will pursue is to employ 
our method for determining the mass scales for the top and W simultaneously in the 
tt di-lepton decay topology. Or, given that the W mass is already quite well-known, 
they could impose this additional constraint in our context and get an excellent t mass 
determination. 

The heart of our technique is the fact that by considering both decay chains in a 
typical LHC event together, a choice for the chain decay masses 7V1 = {my,mx,'m7v} 
(see Fig. ^) in combination with the measured momenta of the 4 visible and measurable 
SM particles emitted in the two chains implies a discrete (sometimes even unique) set 
of three momenta for the two final state A^'s. (One is solving a quartic equation.) 
Conversely, if we have already used our procedure to determine to good precision the 
Ai = {mY,mx,mN} masses, we can invert the process. For each event, we can input 
the known masses and obtain a set of discrete choices for the momenta, pn and pn', of 
the final invisible particles. For each discrete choice, the 4-momenta of all particles in 
the decay chains are then determined. These 4-momenta can then be input to a given 
model (with definite spins for the Y, X, N and definite decay correlations and so forth). 
One can then test the experimental distributions {e.g. of correlation angles, of masses 
constructed from the visible SM particles, and so forth) against predictions obtained 
for the model using a Monte Carlo. Presumably, this will provide strong discrimination 
between different models that have the same already-determined chain decay masses. 
The only question is to what extent the possibility of more than one discrete solution 
for each event will confuse the distributions obtained from the Monte Carlo. 

Conversely, it is clear that determining the spins of all the particles in a chain of 
decays can be difficult without a relatively precise purely-kinematic determination of 
the masses . In particular, we expect that angular correlations and the like (obtained 
from Monte Carlos that assume a particular model including spins) will be strongly 
influenced by the masses. Confusion between two different models with differing spins 
and masses can be anticipated in the absence of an independent purely-kinematical 
determination of the masses. 

Overall, we claim that our techniques provide some powerful new tools for doing 
precision physics at the LHC in an environment where new physics events contain 
invisible particles of unknown mass. We hope the experimental community will pursue 
the approaches we have analyzed. We do not anticipate that fully realistic simulations 
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will lead to significantly larger errors for new particle masses than those we have found, 
but it is clearly important to verify that this is the case. 
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Appendices 



A. Solution Procedure 

To determine whether a given event with the topology of Fig. ^ is consistent with a 
given mass hypothesis, we proceed as follows. We envision the process of 
pp (135) + (246), followed by (135) ^ 5 + (31) and (246) ^ 6 + (42), which in turn 
is followed by (31) — > 3 + 1 and (42) -^4 + 2. (The objects in (. . .) are to be thought 
of as single on-shell particles: in the notation of Fig. |^, (135) = Y, (246) = Y', 
(13) = X, (24) = X', 1 = and 2 = A^'.) We will be assuming input values for m^gg, 
m24g, m^3, m24, and m^, assuming m^gg = ml^Q, ml^ = ml^^ and ml = ml. The 
cross section takes the form 



da 



2s(27r) 



1 



8 




[xiPa + X2Pb - {Pi +P2+P3 + P4 + P5+ Pe)] 



(A.l) 



We first convert 



dxidx2 = -dEtotdpl^ 



'tot ) 



(A.2) 
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introduce on-shell masses for the intermediate particles, and introduce on-shell S 
functions for the invisible particles 1 and 2 to yield 



da 



4(27r) 



1 



8 



J dEtotdpt„tdml^^dml^(idml-^dml^\M\^f{xi)f{x2) 




xd'^Pi6{pl — ml)d'^p2S{pl — ml) 

XS^ [xiPA + X2PB - {Pl + P2 + P3 + P4 + + Pq)] 

x5[{pi + P3 + Pbf - ml^S[iP2 + P4 + PqY - ml^ol 
xS[{p3+Pif - ml^]d[{p2+P4y - 



4(27r) 



1 



8 




xd^Pxb{;p\ — m1)d'^p26{pl — ml) 

x6[{pi + p3 + p^)^ - ml^^]6[{p2 + P4 + Pe)^ - ml^e] 
xS[{p3+Pi)^ - ml^]S[{p2 + P4)^ - ml2] 



(A.3) 



where in the last step we eliminated d^p2 using the 3-momentum conservation part of 
the 5^ function and eliminated Etot using the energy part of the 5^ function. For fixed 
values of the unknown masses, we end up with the 4 unknowns of pf^j and pi, to be 
solved for using the 4 on-shell 5 functions. We will now define 



Assuming no transverse momentum for the y + y = 1 + 2-1-3-1-4-1-5-1-6 system, we 
then have 



Pl ■ P3 = - pIpI - p\pI - pIpI 

P2 ■ P4 = ^2^4 - {plot - Pms ' PDpI ' {'Plis ' PV)PI ' {'Pvis ' Pi)p1 

Pl ■ P5 = E1E5 - pIpI - pIpI ~ pIpI 

P2 • P6 = E2EQ - {pI^ - pI-^ - p'^pI - {-pI-^ - pDpI - (-p^.^ - pDpI . (A.5) 



(Transverse momentum for the Y -\-Y' system can, and must, be included in the 
obvious way. We compute it as the negative of the sum of the observed momenta of 
particles 3, 4, 5 and 6 and the missing momentum.) We next combine the last two 5 
functions and consider the requirement (again, recall that we are assuming some 
input mass values for the intermediate on-shell particle masses) 



Pvis = P3 + P5 + P4 + P6 ■ 



(A.4) 



2p3 • Pl - 2p2 • P4 + A26 = 



Ci^O, 



(A.6) 
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where 

Similarly we combine the 135 and 246 6 functions to obtain 

2pi ■ p5 - 2p2 ■pe + A3b = G2 = 0, (A.8) 

where 

Asb = + ml- mggi + ml^^ - "^42 - "^6 + 2p3 ■ Ps - 2p4 ■ Pe ■ (A. 9) 

Of course, ml and mg are measured experimentally (and are typically small unless 
one is a or Z), and 2p3 ■ p^ and 2p4 ■ are also computable from the experimental 
event. Further, we are assuming input values for m|^, 77142, ''^531 ^^^1 mg42. The above 
is a convenient organization, since /S.2b = and reduces to just momenta dot 
products when the two decay chains are identical. 
We now implement directly the m\i and m\^i 5 functions. 



2pi • p3 = ^li - 2pi ■ p3 = = , (A.IO) 



and 

"^531 - ™ll -^5- 2P3 ■ P5 - 2pi ■ P5 = A531 - 2j9i ■ P5 = G4 = , (A.ll) 

where again 2p3 ■ p^ is determined experimentally and the masses are being input. 
We now solve these 4 equations for the 4 unknowns of p^^^, pi, pf, and pf. We write 
the solutions in the form: 

= C^elEl + 0^.62-^2 + Ccc 
Pi = CyelEl + Cye2-E'2 + Cy 
PI = C^elEl + Cze2E2 + 
plot = CztelEl + Czte2E2 + C^t (A. 12) 

where the c's above are somewhat complicated functions of the masses, energies and 
momenta of the visible particles, 3, 4, 5, and 6. The Jacobian for the variable change 
PiiPiiPiiPtot ~^ Gi,G2,Gs,Gs is easily computed as 

J = 16 [-pspIpIpq + pIpIpIpI + pIpIpIpI - pIpIpIpI 

-pIpIpIpI + pIpIpIpI - pIpIpIpI + pIpIpIp"^ (a. 13) 

It is a function only of observed momenta. For the next stage, we combine the 
expressions of Eq. ( A. 12 ) with 

Pl = -Pl^s-Pl. pI = -pL-pI, P2 = Pit - Plis - pI , (A. 14) 
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and solve the equations for the on-shell 6 functions for pi and p2 (the invisible 
particles) 



= El- {p^f - {p\f - {plf - m\ (A.15) 
^ = El- {plf - {plf - {plf - ml (A. 16) 



for El and i?2- For convenience, we rewrite Eqs. ( [A.15| ) and ( |A.16| ) in the respective 
forms: 

aiiEl + auEiE2 + 022^! + aiEi + 02^2 + a = Fa = (A. 17) 

buEl + buEiE2 + 622^2 + biEi + 62^2 + b = Fb = (A. 18) 

where the Ojj, bij, ai, bi as well as a and b are functions of the c's, mf, m\ and the 
components of Pms- We now take 

F^-^xFb = (A.19) 

Oil 

and solve the resulting linear equation for Ei to obtain 

^ _ Qii & - g bii - a2 bn E2 + an 62 62 - 022 ^ii E^ + an 622 20) 
-aii6i + ai bu + 012 fen -E2 - ctii ^12 E2 



We now substitute this result into Eq. ( [A.17| ) to obtain the final quartic equation for 
E2 of form 

AEl + BEl + CE^ + DE2 + E = 0, (A.21) 

where A, B, C, D and E are functions of the aij, bij, a^, bi, a and b. We then employ 
a standard computer subroutine for obtaining the 4 roots of this quartic equation. 
For typical input visible momenta, some of the roots will be acceptable real solutions 
and some will be imaginary. We retain all real solutions. (The Jacobian for the 
Fa, Fb Ei,E2 transformation is easily computed.) Once real values for Ei and E2 
are obtained these can be substituted into Eq. ( |A.12| ) to determine the 3-vector 
components of pi and the z component of ptot- The components of p2 are then 
obtained by momentum conservation. At this point, the invisible 4-momenta are fully 
determined and could potentially be employed in a model matrix element. 



B. SUSY points 

In this appendix, we give details regarding the SUSY points simulated. 
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B.l Point I 

We input low-scale parameters of 

;x = +300 GeV, tan/3 = 10, (Mi, Ms, M3) = (90, 300, 500) GeV 
~(i,2,3) ^ ~g) ^ ^QQQ Q^y^ ~ (1,2) ^ GeV 

mj^^ = 400 GeV, m[)J^ = 300 GeV, mg^ = m[;% = 1000 GeV (B.l) 

where the m's are the soft slepton and squark masses, and the M's are the gaugino 
masses. L and Q refer to the slepton and squark SU{2)w doublets and E, U and D 
refer to the slepton and squark singlets. Superscripts give the generations. The decay 
chain of interest is 

Ql qxl xl At/2/?, V'R /"X? • (B.2) 



Using the input soft parameters as specified above and SPheno 2.2.3 |^, the 
sparticle masses of relevance for our discussion are (all in GeV): 

m-g ~ 524, mj-^ ~ 438, m^^^c^ ~ 431 

m^o ~ 246.6, m^^ ~ 128.4, m^o ~ 85.3 (B.3) 
For this point, the net cross section available is 

aipp-^ J2 ^L^L+ E ^^^+ E ^^j ~ 2.9x10^ fb, (B.4) 

\ q,q'=u,d,c,s q,q'=u,d,c,s q,q'=u,d,c,s / 

coming from all sources including gg fusion, UlUl fusion, etc. The branching ratios 
relevant to the particular decay chain we examine are 

BR(qL^qxl)r^ 0.27 {q = u,d,c,s) 
BRixl ^ /^H/^^) ~ 0.124 

BRiJil ^ ^^^x\) = 1 ■ (B.5) 

The net effective branching ratio for the double decay chain is 

BR{qLqL ^ 4/ix?X?) ~ (0.27)^ x (0.124)^ ~ 1.12 x lO'^ (B.6) 

for any one qi choice. The effective cross section for the 4/ix'j'Xi final state is then 

(T(4/ix?X?) ~ 2.9 X 10^ fb X 1.12 x 10"^ ~ 32.5 fb . (B.7) 

For an integrated luminosity of L = 90 fb~^, this gives us 2900 A^XiXi events before 
any cuts are applied. After cuts, we are left with about 1900 events. 
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B.2 Point II 

Point II is defined by the following input low-scale SUSY parameters: 

;x = +300 GeV, tan/3 = 10, (Mi, Ms, M3) = (90, 200, 500) GeV 

~(1,2,3) ^ ~g) ^ ^QQQ Q^y^ ~ (1,2) ^ (^^y 

= 400 GeV, m[};5 = 300 GeV, mg^ = m[^% = 1000 GeV . (B.8) 



Using SPheno 2.2.3 ||22[, the relevant chain-decay masses are 

{my = m^o,mx = mji^,mjq = m^o} = {180.8, 147.1, 85.2} GeV . (B.9) 

However, we do not employ the cross sections and branching ratios predicted by these 
parameters. Instead, we assume 2000 available experimental points after cuts, close to 
the 1900 left after cuts in the case of Point I. This allows us to see how errors change 
in the case where my — mx is much smaller than in the case of Point I. 

B.3 Point III 



The masses used in this case are obtained from PYTHIA 1.0.8 using the low-scale 
parameters 

^ = +3000 GeV, tan/3 = 10, (Mi,M2,M3) = (0.2,200,500) GeV 
= ^g) = 1000 GeV, m^^'^^ = 100 GeV 
= 400 GeV, mj);5 = 300 GeV, mg^ = m[;% = 1000 GeV (B.IO) 

yielding 

{my = m^o, mx = mji^,mx = "^x?} ~ {199.4, 100.0, 0.1} GeV . (B.ll) 

Again, we assume 2000 available experimental points after cuts, close to the 1900 
events after cuts obtained for Point I. 

B.4 Point IV: SPSla 

For the SPSla point, we use the GUT-scale mSUGRA inputs of 

mi/2 = 250 GeV, mo = 100 GeV, Aq = -100 GeV, tan/5 = 10, fi > 0. (B.12) 
From IS A JET 7.75, the spectrum is calculated as 

rUg ~ 608 GeV, m^^ ~^ ~ 571 GeV, ms^,^^ ~ 565 GeV 



m^o ~ 180.3 GeV, m^^ ~ 142.5 GeV, m^o ~ 97.4 GeV rrir, ~ 134.7 GeV 



(B.13) 
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The total effective cross-section including all channels for X2X2 production is about 
1 pb. The relevant branching ratios are: 

5^(X^ ^ ~ 0.063, Si?(/I±^/.±x?) = l- (B.14) 

For L = 300 fb^^, this gives 1200 events before any cuts. After cuts, we are left with 
about 425 events. Errors for the masses given in the text are based upon the latter. 
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